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Abstract 

The sampled negative energy density seen by inertial observers, in arbitrary quantum states is 
limited by quantum inequalities, which take the form of an inverse relation between the magnitude 
and duration of the negative energy. The quantum inequalities severely limit the utilization of 
negative energy to produce gross macroscopic effects, such as violations of the second law of ther- 
modynamics. The restrictions on the sampled energy density along the worldlines of accelerated 
observers are much weaker than for inertial observers. Here we will illustrate this with several 
explicit examples. We consider the worldline of a particle undergoing sinusoidal motion in space 
in the presence of a single mode squeezed vacuum state of the electromagnetic field. We show 
that it is possible for the integrated energy density along such a worldline to become arbitrarily 
negative at a constant average rate. Thus the averaged weak energy condition is violated in these 
examples. This can be the case even when the particle moves at non-relativistic speeds. We use 
the Raychaudhuri equation to show that there can be net defocussing of a congruence of these 
accelerated worldlines. This defocussing is an operational signature of the negative integrated en- 
ergy density. These results in no way invalidate nor undermine either the validity or utility of the 
quantum inequalities for inertial observers. In particular, they do not change previous constraints 
on the production of macroscopic effects with negative energy, e.g., the maintenance of traversable 
wormholes. 
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I. INTRODUCTION 



It is well known that quantum field theory allows for the existence of negative energy 
density, which constitute local violations of the weak energy condition. For a recent review, 
see Ref. pQ. Negative energy density can arise either from boundaries, as in the Casimir 
effect, from background spacetime curvature, or from selected quantum states in Minkowski 
spacetime. The last possibility will be the focus of the present paper. It is possible to create 
states, such as a squeezed vacuum state of the quantized electromagnetic field, in which the 
energy density at a given spacetime point is arbitrarily negative. However, the duration 
of the negative energy is strongly constrained by quantum inequalities [2HH]. These are 
restrictions on a time averaged energy density measured by an observer. (Time averaging 
is essential, as there is no analogous restriction on spatial averages [TO]-) Let us consider 
the case of inertial observers in Minkowski spacetime, with four velocity u** '. If (T^) is 
the expectation value of the normal ordered stress tensor operator in an arbitrary quantum 
state, then quantum inequalities take the form 



Here r is the observer's proper time, /(r) is a sampling function with characteristic width 
ro, and d is the number of spacetime dimensions. The dimensionless constant Co depends 
upon the form of the sampling function, and is typically small compared to unity. In the 
limit r — > oo, Eq. ([TJ becomes the averaged weak energy condition 



which states that the integrated energy density along an inertial worldline is non-negative. 
The essence of a quantum inequality is that there is an inverse relation between the mag- 
nitude and duration of negative energy density. These relations place strong constraints on 
the effects of negative energy for violating the second law of thermodynamics [2], and for 
maintaining traversable wormholes [TT] or warpdrive spacetimes [12J. 

A more general quantum inequality for arbitrary worldlines has been proven by Few- 
ster [13J. However, this inequality is often very difficult to evaluate explicitly and can be 
very weak. There are some known examples where the integrated energy density along a 
non-inertial world line can be arbitrarily negative. One example comes from the Fulling- 
Davies moving mirror model in two spacetime dimensions [HI [15] . A mirror with increasing 




(1) 




(2) 
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proper acceleration to the right can emit a steady flux of negative energy to the right. An 
inertial observer could only see this negative energy for a finite time before being hit by the 
mirror, and the integrated energy density seen will be consistent with Eq. ([I]). However, an 
accelerated observer who stays ahead of the mirror can see an arbitrary amount of negative 
energy This example suffers from two unrealistic features: it can only be formulated in two 
spacetime dimensions, and it requires an observer with ever increasing proper acceleration. 

A second example was provided by Fewster and Pfenning [16J, who analyzed the case of a 
uniformly accelerating observer in the Rindler vacuum state. This state has negative energy 
everywhere within the Rindler wedge. An observer with constant acceleration can also see 
an arbitrary amount of negative energy. However, the constant acceleration requires the 
observer to move arbitrarily close to the speed of light and hence have an unlimited source 
of energy. It is also not clear whether the Rindler vacuum is a physically realizable state. 

The main purpose of this paper is to construct some more realistic examples of accelerated 
motion in which the observer can have arbitrarily negative integrated energy density. We 
will consider observers who undergo sinusoidal motion in the presence of a squeezed vacuum 
state of the quantum electromagnetic field. We find that even in the case of non-relativistic 
motion, it is possible for the integrated energy density in such an observer's frame to grow 
negatively at a constant rate in time. In Sect. [TTJ we consider a squeezed vacuum state for 
a single plane wave mode, and motions both perpendicular and parallel to the direction 



of propagation of the wave. In Sect. Ill, we repeat the analysis for the lowest mode in a 
resonant cavity in a squeezed vacuum state. In Sect. [TV] we address a possible physical 
effect of accumulating negative energy density, in the form of defocussing of a congruence of 
accelerated worldlines. Our results are summarized and discussed in Sect. |Vj In particular, 
we argue that the results in this paper neither contradict, nor diminish the utility of, the 
usual quantum inequalities proven for inertial observers. 

Throughout this paper, units in which h — c—1 will be used. Electromagnetic quantities 
are in Lorentz-Heaviside units. 



II. OSCILLATIONS THROUGH A PLANE WAVE 

Let us first evaluate the stress tensor components for a single mode plane wave in a 
squeezed vacuum state of the electromagnetic field. The electromagnetic stress tensor is 
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given in terms of the field strength tensor as 



T a p = F ap Fp p — -g a /3 F^F^ . (3) 



Its spatial components are 



Tjl = _ EjEl _ BjBl + 1^ (E 2 + B 2 ) , (4) 

the energy density is 

T tt = -(E 2 + B 2 ), (5) 
and the energy flux in the i-direction is 

T u = (E x B)* (6) 

Write the electric and magnetic field operators in terms of photon creation operators a kA 



and annihilation operators a kA as 



E = ^(a kA ^ kA + 4 A r kA ), (7) 

k,A 

and 

B = £(a kA £ kA + 4 A £* kA ). (8) 

k,A 

Assume that the excited mode is a plane wave propagating in the z-direction, with polar- 
ization in the x-direction. Then its mode functions take the form £ kA = x£, and £> kA = y B, 
where 

S = B = J^-e^'- t K (9) 
V £ v 

Here V is the quantization volume and Q = |k| is the angular frequency of the wave. 
Quadratic operators are assumed to be normal ordered with respect to the Minkowski vac- 
uum state, so 

(EjEt) = 5 jx 5 lx (E 2 ) = S 2 (a 2 ) + (S*) 2 ((a^) 2 ) + 2\S\ 2 ^a) . (10) 

where a is the annihilation operator for the excited mode. Similarly, 

(B ] B l } = 5 jy 5 ly (E 2 ). (11) 

The quantum state is taken to be a single mode in which case 

(E 2 ) = (B 2 ) = 2fte[sinh 2 r \S\ 2 - £ 2 sinhr coshr e* 5 )] 

= ^ sinhr jsinhr — coshr cos [2Q(z — t) + 5} | , (12) 



where r is the "squeeze parameter" and 5 is a phase parameter. The nonzero components 
of the stress tensor are given by 

(T a ) = (T zz ) = (T tz ) = (E 2 ) . (13) 



We see from Eqs. (12) and (13) that the energy density can be periodically negative in the 
lab (i.e., inertial) observer's frame, but the positive energy density always outweighs the 
negative energy density, in accordance with the quantum inequalities. 

The energy density in the inertial frame has its minimum (most negative) value when the 



cosine term in Eq. (12) is one, so 



(T") > sinhr (coshr - sinhr) = (1 - e 2r ) > ■ (14) 

Thus the maximally negative energy density is bounded below, and occurs for large r. 
However, in this limit the maximally positive energy density is unbounded and grows as e 2r . 
In the opposite limit, where r < 1, the energy density is approximately oscillatory 

(T tt ) « -r ^ cos[2ft(z - t) + 5} + r 2 ^ + 0(r 3 ) . (15) 

However, there is also a positive non-oscillatory term of order r 2 . 



A. Perpendicular Motion 

Now consider a non-geodesic observer who moves on a path which is perpendicular to the 
direction of propagation of the wave. Let this path be defined by 

(1 T 

v x {t) = ^ = Asm{ut), (16) 

where \A\ < 1, and v y = v z = 0, and u is the angular oscillation frequency of the observer's 
motion, and where we have chosen z = 0. Then 

l = ^= 2 = - r 1 9 > (17) 
and the observer's four-velocity (as measured in the lab frame) is 

= 7(1, v x ,0,0), (18) 

where u* = 7 = dt/dr. 



The integrated energy density along the accelerated observer's worldline is 

/ = j(T^u u )dr, (19) 

where the integrand is 

(T^u^)dr = 7 2 (T u )dr = (T a ) . dt . (20) 

A - A 2 sm 2 (ut) 



Here we used the facts that (T tx ) = (T xx ) = and 7 2 <ir = jdt. If we expand to first order 
in r, the result is 

(T^u,u v )dr - - r J^£BzM . (21) 
Vyjl - AHm 2 {ut) 

The numerator of this expression describes the fact that, for small squeeze parameter, the 
inertial frame stress tensor components are nearly sinusoidal. The denominator describes the 
effect of going to the non-inertial frame. If we can arrange that the 7 factor has its maximum 
value when the numerator is negative, then accelerated observer will see net negative energy. 
This situation occurs when u> = Q and when 5 = tt, as illustrated in Fig. [T] We will make 
this choice throughout the remainder of this subsection. 
In this case, the integrated energy density becomes 

rQ f dt / ^ s 

cos(2fit). (22) 



V J yfl - A 2 sin 2 {Qt) 
If we perform the integration on t and multiply by the quantization volume, we get 

l v ^ r [2 E(Qt, A 2 ) + {A 2 - 2) F(Q t, A 2 ) ^ ^ 

where F(Qt,A 2 ) and E(Qt,A 2 ) are elliptic integrals of the first and second kind, respec- 
tively. 

As a specific example, let us plot IV for r = 0.01, A = 0.9, and in units where Q = 1. 
Since, strictly speaking, the energy density is inversely proportional to V, we want to make 
a graph of I V as a function of r, i.e., a graph of the integrated energy density, multiplied by 
the quantization volume, seen by the accelerated observer as a function of his proper time. 
The relation between r and t is r = / dt/7, which is 

E{Qt,A 2 ^ 



n 



(24) 



If we plot Eq. (23) against Eq. (24) for our chosen parameters, we get Fig. 2 



FIG. 1: The figure illustrates that maximum negative energy density is obtained when we set uj = Q 



and 5 = ir. The dotted line represents the Lorentz factor in Eq. (21), [1 — A 2 s'm 2 (Q t)] -1 / 2 , while 
the solid line represents the cosine term, — cos(2 Q t — 5) = cos(2 Vt t) both graphed as functions of 
time. Here we have chosen A = 0.2 and Q = 2. (The figures have been appropriately scaled to 
allow easier visualization.) 

Now let us examine our expression for / in the A <C 1 limit. If we expand the Lorentz 
factor to second-order in A, we obtain 

1 i 

l2 • 2/ojA i sn/ /l4 



1 - A 2 sin 2 (ttt) 



2/o^ 2 



1 + -A 2 sm 2 {Q t) + 0{A A ) . (25) 



In this limit, the difference between dt and dr will be 0(A 2 ). If we use Eq. (25) in Eq. (22) 
to calculate J, we find: 

rA 2 VLT rsin(2ftT) r A\m(2VLT) r A 2 sin(40T) , , 

I ps 1 H i i . (26) 

8V 2V SV 32V y 1 

The sinusoidal terms will eventually be dominated by the linear term, but this can take 
many cycles, so we keep the O(A ) sinusoidal term, but drop the 0(A 2 ) sinusoidal terms. 
Therefore, our two leading order terms are 

rA 2 VLT rsin(2fiT) , . 

I « + v . (27) 

8V 2V K ! 
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FIG. 2: The integrated energy density multiplied by the quantization volume, IV, seen by an 
accelerated observer who is moving perpendicularly to the direction of wave propagation, is plotted 
as a function of his proper time, r, for the parameters r = 0.01, A = 0.9, and in units where Q. = 1. 

Here A 2 <C 1, so the oscillating term is larger in mag nitude until T > A/(A 2 VL). After 
this, the linear term dominates. However, we should recall that there is positive r 2 term in 



Eq. (15). This term will give a contribution to / of VLi T/V, and is negligible only if we 



require that 

8r<A 2 . (28) 
Nonetheless, accumulating negative energy density, can occur for arbitrarily small velocities. 



For any A ^ 0, we can find a value of r which satisfies Eq. (28). Then eventually the first 



term in Eq. (27) will dominate 



B. Parallel Motion 



We now consider the case of the accelerating observer moving parallel to the direction 
of the propagation of the wave. In the lab frame, we have (T tt ) = (T zz ) = (T tz ). The 
accelerated observer's three-velocity and position, respectively, are 

dz 



v*(t) 



dt 



Asm.{uit) , 



(29) 



z(t) 



A 



cos (a; t) 



and so 



u" = 7(1,0,0,*;) 



where it* = 7 = dt/dr. Therefore, we have that 



(T t ' t ')= 1 2 (l-2v + v 2 ) (T tt ) 



1-v 



(T tt ) 



(30) 



(31) 



(32) 



where (1 — v)/(l + v) is a linear Doppler shift factor (as opposed to the transverse Doppler 
factor in the perpendicular case). As a result, 



(T^)u.u v dr 



1-v" 
1 + v 



(T tt ) dr 



X3/2 



vTT 



(T**> dt 



(33) 



since dr = 7 1 dt = yl — t> 2 <it. 

Here the observer is moving in the direction of wave propagation, so we can no longer 



set z = 0. Now the energy density in the inertial frame is given by Eq. (12), with z 
— (A/cj) cos(a;t), so 



(T tt ) 



o 
v 



sinh r < sinh r — cosh r cos 



/2AQ\ 
\ uj 



\ cos(wt) + 2Qt - 8 J . 



(34) 



In this case, we find accumulating negative energy density for 00 = 2 0, and 5 = —ir/2. The 



integral in Eq. (33) can be done analytically for small A, as will be discussed below, but 



for more general A, it can only be performed numerically. As an example, let us choose the 
case where r = 0.01, 5 = —n/2, A = 0.9, and to = 2, in units where = 1. In Fig. |3j we 



graph IV against the observer's proper time, which will again be given by Eq. (24). 



Now we wish to consider the non-relativistic limit, and work to first order in v and hence 
in A. To this order, dr f» dt, so 

[(T^dt, (35) 
where the energy density in the accelerating frame is 



(T*'*') sa (1 - 2v) (T u ) . 



(36) 



If we expand Eq. (34) to first order in A, the result is 



(T tt ) P3 — sinhr <! sinh r - cosh r 



20 

cos(20t - 5) A cos(ut) sin(20t - 5) 

UJ 



(37) 
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FIG. 3: The integrated energy density multiplied by the quantization volume, IV, seen by an 
accelerated observer who is moving parallel to the direction of wave propagation, is plotted as a 
function of his proper time, r, for the parameters r = 0.01, A = 0.9, and u = 2, in units where 
fi = 1. Note that the accumulated negative energy density grows much faster than in the case of 
the perpendicularly moving observer, due to the linear Doppler shift term in the energy density. 
We also see extra structure in this curve as well, because the expression for the energy density is 
more complicated than in the perpendicular motion case. 



where we have used the fact that 

'2VLA 



cos 



sm 



CO 



cos{ojt) 
cos(ut) 



1 + 0(A 2 
(2VlA 



uj 



cos(wi) + 0(A 2 



(38) 



Next we evaluate the energy density in the accelerating frame to first order in A and set 
cu = 2Q to find 



(T ) 



n 



v 

coshr 



sinhr {sinhr [1 — 2Asm(2Qt)} 



3 1 

cos(2^t - S) Asm(4ttt -5) A sin<5 

2 2 



(39) 



This expression reveals that we can have growing negative energy density if 5 = —it/ 2 and 
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r -C 1. In this case, we may write 

O r _ n \i 

(40) 



r 2 -r(-A- sin(2fit) 



where order A oscillatory terms have been dropped. If 



r<.^A, (41) 



which is the analog of Eq. (28), the integrated energy density grows negatively as 

rAQT r . . r AVtT . , 

+ — [l-cos(2fiT)] — . (42) 



2V 2V l v /J 2V 

The latter asymptotic form holds for 

T>> ^- < 43 > 

In the parallel motion case, the rate of growth of the negative integrated energy density is 
first order in A, as compared to second order in the perpendicular motion case treated in 
the previous subsection. This is due to the fact that in the parallel case, there is a linear 
Doppler shift, whereas in the perpendicular case the Doppler shift is transverse. 



III. OSCILLATIONS IN A CAVITY 



A. The Perpendicular Case 



We now consider the case of a particle oscillating in a closed cavity with dimensions a, b, 
and d aligned along the respectively, where b < a < d. The modes in this cavity 

were discussed in Ref. [17J. With the condition that b < a < d, the lowest frequency mode 
is the TE mode with p — I — 1, m — 0, where the frequency of the mode is given by 



n = -V? + ?- (44) 



and the non-zero components of the electric and magnetic fields are 

£ x £ z , 

£ v = — C sm(-x) sin (~z) e~ lil1 



ix \a ' \d 

B x = i — C sin \ —x) cos \ —z] e 



i 



By = 



B z = -iC cos ( n x) sm( 7r -z)e- in \ (45) 
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where the electric field is taken to be polarized in the y-direction. This mode is independent 
of y. Here C is a real normalization constant, given by 

2VL 



C z = 



(46) 



abd(l + a 2 /d 2 ) ' 

For the case where only a single mode j is excited, the normal ordered expectation values 
of the squared fields are 



and 



where 



(E 2 ) = 2(a t a) \Sj\ 2 + 2Re^(a 2 } £ 2 ^ 



(B 2 ) = 2(a^ a) \Bj\ 2 + 2Re((a 2 ) B 2 



(47) 



(48) 



2„2 



c2 _ c 2 _ " " 



C 2 sin 2 (-:r) sin 2 (J^z^j e 



-2iM 



-c 2 



2/ 71 " \ • 2/ 71 " \ ° • 2/ 71 " \ 2/ 

cos —x sin ( ) H — — sin — x cos 
v a ^ v a 7 a 2 v a ^ v 



'7T 



-2iQi 



(49) 



In this case, we can have the particle moving in the y-direction, and located in the center 
of the cavity in the other directions, so that 

a d 



x = -, z = - 
2 2 



This considerably simplifies the mode functions, leading to 

fla 



(50) 



B x = B z = 0, 



Ce 



7Y 



(51) 



The only non-zero components of the stress tensor which we will need are (T tt ) and (T yy ), 
which become 



(T* 



■(T 



yy\ 



sinhr 



sinhr \S y \ 2 — cosh r Re {e lS S 2 ^ 



= N sinhr [sinhr — coshr cos(2f2t — 5)] , 



where 



7T 



Because the direction of oscillation of the particle is in the y-direction, 



«" = 7 (1, 0,^,0), 



(52) 



(53) 



(54) 
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where v y = A smut. The integrand of I is 



(T. u u^)dr = iilkl±|£l£k^ dt = J I - A 2 s\n 2 (ut) (T tt ) dt . (55) 
1 - A 2 sm 2 (cut) 



Let 



U= Jl-A 2 sm 2 (ut)(T tt ). (56) 



We will assume that A -C 1, and expand to second order in A. Therefore, if we use Eq. (52) 



in Eq. (55), and set u = f2, we have that 



U « N (l - - A 2 sin 2 (f2t)) [sinh 2 r — sinhr coshrcos(2fii — 5)] . (57) 



If we now expand the right-hand side of Eq. (57) to second order in r, set 6 = 0, integrate 
from to T, and drop oscillatory terms in A 2 , we obtain 



I~N[r 2 T-\A 2 rT-r ^ T) Y (58) 



where we have also dropped a higher order A 2 r 2 term. 

The positive first term is negligible compared to the second when 

r < I A 2 . (59) 
8 

The middle negative linear growing term will dominate the sinusoidal term when 

T > At < 60 > 

In this case, the restrictions on r and T are the same as those for perpendicular motion in 
the plane wave case. In these limits, we therefore have negative energy density which grows 
linearly as 

/ « — - NA 2 rT. (61) 



If we use Eqs. (46) and (53), we can write the previous equation as 

rA 2 QT 



(62) 



AV 

where V = abd is the volume of the cavity. Compare this result with the first term in 



Eq. (27), the corresponding rate for perpendicular motion in a plane wave mode. If we 
identify the cavity volume in the former with the quantization volume in the latter, then 
they differ only by a factor of two. 
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B. The Parallel Case 



In this subsection, we will consider the case of a particle oscillating in a cavity along the 
z-axis, in the limit where A -C 1, and work to first order in A. We take 



v 7 = A sin cut 



(63) 



and 



Z = Z 



A cos ut 



10 



+ 0(A 2 ), 



(64) 



where z = zq corresponds to the equilibrium position of the particle, and the last term 
corresponds to relativistic corrections. We will choose the x-position of the particle to be 



x 



(65) 



The energy density in the particle's frame is 



(T*'*') = 7 2 ((T tt - 2v z T tz + v 2 T zz )) « (T tt ) - 2v z (T tz ) 



(66) 



Here 



(T tz ) = -(E y B x ) 



(67) 



We need to calculate (T a ) and (T tz ) using the mode functions in Eq. (45), and then expand 
the result to second order in r. The result for (T* * ) may be written as 



(r* * } « C 2 {(Fi + 2v z F 3 )r 2 - r [F 2 cos(2fit - 5) + 2v z F 3 sin(2fit - 5)}} , (68) 



where 



and 



F 1 = F 1 (z) 
F 2 = F 2 (z) 



fl 2 a 2 

n 2 a 2 

7T 2 



Fn(z) 



TCZ 



TTZ 



siir m + ^ cos u 



7rz \ a 



sin 



TXZ 

T 



2tt2 



Qa 2 
lid \ d 



and where C 2 is once again given by Eq. ( 46 ) . 



The integrated energy density may be written as 



Injdt (T*'*')~C 2 (r 2 h-rl 2 



(69) 
(70) 

(71) 



(72) 
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where 



and 



h = I* dt (F, + 2v z F 3 
Jo 



(73) 



J 2 = ( T dt \F 2 cos(2ftt - 5) + 2v z F 3 sin(2fit - 5)] . (74) 
Jo 

As in the case of parallel motion in the plane wave case, with the appropriate choices for 
oo, f2, r and 5, we expect to get a linearly growing negative term, a term which is first order 
in r and sinusoidal in time, and a positive r 2 term. The first and second of these terms will 
arise from F 2 and F%, while the third term will arise from F±. We also expect that we will 
find a non-trivial effect in first order in A. 



Let us first examine the terms involving F 3 in Eq. (68). These terms both involve the 



product v z F 3 , and are hence already of order A. Thus we may use Eq. (64) to write 

Qa 2 



F«(z) 



F 3 (z ) 



n d 



(75) 



where we have set zq = d/4. (As it turns out, the linearly growing term we want will come 
from the F 3 term in 1%, so we cannot choose Zq = d/2.) 

A similar situation applies to Fi, which contributes only to an order r 2 term. This is a 
positive, growing term which we need only to zeroth order in A. For this purpose, we may 
evaluate F\ at z — Zq\ 

F 1 (z)^F 1 (z ) 
Thus, for estimating the order r 2 term, we may use 

h~F 1 T, 



fl 2 a 2 
2tt 2 



1 + 



d 2 



(76) 



(77) 



where F\ has the value in Eq. (76). 



The negatively growing term comes from J 2 , which involves F 2 , so we need to expand the 



latter to first order in A, using Eq. (64), as 

n 2 a 2 



Foz) 



2tt 2 

n 2 a 2 



1 - 



2tt 2 



d 2 / 
a 2 / 
d 2 



a 

~d 2 i 
„2 S 



■2tiz\ 



1+ ^J (^) cos( ^ 



(78) 



The I 2 term will be maximally negative when u = 2VL and 5 = 0. In this short 
calculation yields 



, ^a 2 / 
And \ 



a 
d 2 



T 



Qa 2 
Air 



1 



d 2 



sin(2fiT) 



(79) 
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where oscillatory, order A terms have been dropped. Note that 3 — a 2 /d 2 > because a < d. 
Therefore, the integrated energy density becomes, 



C 2 



-rAT 



Qa 2 
And 



3 - 



d 2 



— r 



Qa 2 
Jtt 2 



1 



a 
d 2 



n 2 a 2 

~2^ 



sin(2fiT) + r 2 T 



1 + 



d 2 



(80) 



We see that the negative linearly growing term will dominate the sinusoidal term when 

m d ( d 2 — a 2 
T > 



ttA \3d 2 -a 2 



and the positive, order r 2 term, when 



2nd ( d 2 + a 2 



r < A. 



7T 



(81) 



(82) 



3d 2 -a 2 , 

In this case, we find that the integrated energy density in the particle's frame grows nega- 
tively as 



rAQ 
^2V 



T 



Qa 2 3d 2 -a' 



nd \ d 2 + a 2 



(83) 



where we have used the definition of C 2 and the fact that V = abd is the volume of the 



cavity. Compare this result with Eq. (42), the corresponding rate for parallel motion in 
a plane wave mode. If we identify the cavity volume in the former with the quantization 
volume in the latter, then they differ only by the factor in the square brackets. If a and d 
are of the same order of magnitude, then Eq. (44) tells us that Q ~ 0(l/a) ~ 0(l/d), and 
this factor is of order unity. 



IV. EFFECTS OF THE NEGATIVE ENERGY ON FOCUSSING 

In this section, we will treat one possible effect of the accumulating negative energy along a 
particle's worldline. It is well-known that the attractive character of gravity, with ordinary 
matter as a source, leads to focussing of null and timelike geodesies. One expects that 
negative energy densities might have the opposite effect, and produce defocussing through 
repulsive gravitational effects. 



A. Raychaudhuri Equation 



The effect of gravity on a congruence of timelike worldlines is described by the Raychaud- 
huri equation. In our case, we allow the worldlines to be non-geodesics, so the equation takes 
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the form HE 



d8 1 

= — = -R a pu a u? + 2Lo afS uj a P - 2a aP (T aP - -6 2 + V p aP . (84) 

GST 3 

Here u a and = u a V 'aU 13 are the 4-velocity and 4-acceleration of the congruence, and 
a a R and u a s are the shear and vorticity tensors. Also, 9 = V a u a is the expansion, and 



R a p is the Ricci tensor. The last term in Eq. (84) is the acceleration term, which vanishes 
for geodesies. We will assume a hypersurface orthogonal congruence, in which case the 
vorticity tensor vanishes, uj a p = 0. In addition, we assume that the shear and expansion 
are sufficiently small, that the terms quadratic in those quantities may be neglected. In this 
case, the Raychaudhuri equation becomes 

6 W -J^uV + ^ac , (85) 

where (9 ac = V 'pa 13 is the acceleration term, and the Ricci tensor term describes the effects 
of gravity. 

Next we assume that an electromagnetic field is both the cause of the acceleration and 
the sole source of the gravitational field. Particles with rest mass m and electric charge q 
obey the equation of motion 

ap = ^-Fp p u p , (86) 
m 

where the field strength tensor, Fp p , is assumed to obey the source free equation 

V Q F aP = . (87) 

We can now write the acceleration term as 

e\ c = ^-F a \V a u p ). (88) 
m 

The covariant derivative of the 4-velocity may be expressed as [19] 

V a up = ap a + \^9(g a p + u a up) - a p u a , (89) 

when up a = 0. However, all terms on right hand side of this expression, except for the last, 
are symmetric tensors which vanish when contracted into the antisymmetric field strength 
tensor. Thus we obtain 

0\ c = -apa? . (90) 
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The electromagnetic stress tensor, given in Eq. ^ is tracefree, so the Einstein equations 
become 

R Qf} = 8ne p 2 T^, (91) 
where £ p is the Planck length, and Newton's constant is G = £ p 2 , in units where h = c = 1. 



We may write 



1 771 ^ 1 

T a/J uV = {u a F ap )u^F^ + -F, U F^ = — a p aP + -F, U F^ , (92) 



where we have used Eq. (86). We may use this expression to evaluate the Ricci tensor term 



in the Raychaudhuri equation and write 



0«-fl + 87r^^aX - 2tt£ p 2 F^F^ . (93) 
B. Fields Producing Acceleration 

In previous sections, we assumed a prescribed sinusoidal motion, but did not explicitly 
give the electromagnetic fields which would produce this motion. Here we will concentrate 



on the case of motion parallel to a plane wave mode, which was treated in Sec. II B In 



particular, we consider the case of non-relativistic motion along the ^-direction, as described 



by Eq. (29), with A < 1. This motion can approximately be produced by a plane wave 
with polarization in the z-direction. Here we will consider a classical electromagnetic wave 
propagating in the ^-direction, with electric field E = E c z, and magnetic field B = E c ~k, 
where 

,_, lo Am , . . 
E c = - cosu(t - y) . (94) 

To order A, only the electric field determines the motion of the particle, with the magnetic 
force contributing in order A 2 . Because the motion of the particle is only in the z-direction, 
we may set y — 0. In this case, the energy density of the classical wave, in the laboratory 
frame, is 

T? = E 2 C = (^^J cosV). (95) 

In addition to this classical field, the particle is also subjected to the quantum fields asso- 
ciated with the squeezed vacuum state mode. These fields potentially produce a fluctuating 
force on the particle, which we wish to include. Let E q and B g be the terms in Eqs. ^ and 
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QSJ) , respectively, which refer to the mode in a squeezed vacuum state. That is, 

E, = x(£ a + £*a f ), (96) 

and 

B q = y(S a + £*if), (97) 

where £ is defined in Eq. ([9]). We will treat the velocity of the particle due to the quantum 
electric field as an operator in the photon state space, v 9 = v q x, where v q will be evaluated 
explicitly below. 

There is a third effect which we will not include explicitly. This the effect of the emitted 
radiation by the particle. There will be an average radiation reaction force which will slightly 
change the trajectory of the particle for a given classical field. However, this is normally 
very small and will be neglected. There will also be a shot noise effect, an uncertainty in the 
particle's momentum due to the statistical uncertainty in the number of photons emitted. 
This effect depends primarily on the classical field driving the average motion and not upon 
the quantum electric field. Hence it, and the radiation reaction force, would cancel in any 
experiment which compares particle motion with and without the quantum electric field. 
In addition, this momentum uncertainty grows as the square root of the mean number of 
photons radiated, and hence as the square root of time. Here we are interested in effects 
which grow linearly in time. 

We will now compute the components of the acceleration four-vector in the lab rest frame, 
taking account of both the classical and quantum parts of the electromagnetic field and of 
the particle's four-velocity. The acceleration four-vector satisfies 

a P = <L pp<* Ua . (98) 
m 

In the non-relativisitic limit, the four- velocity is 

u a = {l,v q ,0,v c ), (99) 
where v c = A sin(o;t) . The non-zero components of the field strength tensor are 

F tz = pyz = ^ ptx = pzx = p^ ^ ( 10 q) 

and those obtained by antisymmetry of F pa . The components of a p become 

a* = ^(E q v q + E c v c ) 



m 
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a x = Q-E q (l-v c ) 
m 

a y = — E c v c 
m 

a z = ^(E qVq + E c ). (101) 
m 

We can now form the scalar a p a p , expand it to first order in the velocities, dropping v 2 , 
v 2 and v c v q terms, and take its expectation value in the squeezed vacuum state. The result 
is 

a 2 r i 
(a p a p )^\ (E 2 )(l-2v c ) + E 2 c +2E c (E q v q ) . (102) 

m z l « J 

Let us examine each term on the right-hand side of this expression. The classical energy 
density, which is the same to first order in velocity in the lab frame and in the particle rest 
frame, is just E 2 . Because the classical wave is propagating in the y-direction, and all the 
motion is in the x and z directions, this is the perpendicular motion case, with respect to 
the classical wave. Thus, to order v, T l c 1 w T", since 

T t't> w T tt + o(v 2 ). The expectation 
value of the quantum energy density in the lab frame is (E 2 ), and is given explicitly by 



Eq. (12). This quantity in the particle rest frame is {E 2 )(l — 2v c ). The final term is the 
contribution of the velocity fluctuations to the acceleration. 

For both the classical and quantum electromagnetic fields, we have assumed plane waves, 



for which E 2 = B 2 , and hence F llv F pl ' = 0. Thus we may drop the last term in Eq. (93), 



and write mean rate of change of the expansion as 



Cm 2 " 

,2 



C. Velocity Fluctuations and Defocussing 



(0)«- 1 + 8tt^— (a/). (103) 



The fluctuating part of the velocity, v q , is determined by Eq. (101): 



*fc = „ = ^(l-*), (104) 

where the term proportional to v c on the right hand side is due to the magnetic force 
produced by B y . Note that time derivative here is a total derivative, and we need to account 
for both the explicit time dependence and the implicit dependence through z(t): 

dv 1 _dv q dv^ 

dt~ dt + dz Vc > (iU5j 
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recalling that v c = dz/dt. The solution to Eq. (104) becomes 



iq 



{Sa-S* a ] ) 



(106) 



q m n 

where we have used Eqs. (|9| and (96). Note that the effects of the magnetic force and of 
the implicit time dependence cancel one another. 

We may compute (E q v q ) = (: E q v q :) in a squeezed vacuum state to find 

Q 



(E q v q ) 



mV 



sinhr coshr sin[2fi(,z — t) + 8] . 



Here we used 



(a 2 ) = ((a 1 ") 2 )* = -e ld sinhr coshr 



(107) 



(108) 



in the squeezed vacuum state. We may use Eq. (J94J) with y — 0, and set 8 = —tt/2 to write 

(109) 



2E C (E q v q ) 



2uA 



V 



sinhr coshr cos(u;t) cos[2fi(z — t)] . 



We will work only to first order in A, which means that we can ignore the ^-dependence 



(see Eqs. (30) and (38)) in the above expression, which will contribute in order A 2 . When 
we set u = 2Q, and drop oscillatory terms, then we have 

. . 2QA 
2E C \E q v q ) ~ ti— sinhr coshr . 



V 



In the small r limit, this becomes 



2E C {E q V q ) 



2VLA 
V 



(110) 



(111) 



which is to be compared with the same limit for the squeezed state energy density in the 
accelerated frame, 

O A 

(112) 



(E 2 q )(l-2v c )^-^r. 



2V 



The latter quantity is just the order r, non-oscillatory term in Eq. (40). We see that both 
terms have the same form and same sign, and both contribute to defocussing, although the 
effect of the quantum velocity fluctuations is four times that of the negative energy density 
in this limit. 

If we combine these terms, as well as the time average of the classical energy density, 



Eq. (|95|), evaluated at oj = 2fl, then Eq. (|102|) for the mean squared acceleration becomes 

(113) 



(a p a p ) w 2 n 2 A 2 - 



5q 2 ttAr 



AVm 2 
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The positive term is the focussing effect of the classical energy density, and the negative 
term is the combined defocussing effect of the negative energy density and the velocity 
fluctuations. These two terms depend upon different combinations of parameters, and it 
seems possible to arrange for the defocussing effect to dominate. Note that the gravitational 



effect, from the Ricci tensor, is oc £ p 2 in Eq. ( f 03 ). The part without £ p z is a pure acceleration 



2 • 



effect from the acceleration term. However, both effects have the same functional form here. 

V. SUMMARY AND DISCUSSION 

The key result of this paper is that an accelerated observer undergoing sinusoidal motion 
in space can observe an average constant negative energy density, so the integrated energy 
density grows negatively in time in this observer's frame. This is contrast to an inertial 
observer, in whose frame the energy density is more constrained by quantum inequalities. 
We considered a squeezed vacuum state for both a plane wave and a standing wave in 
a cavity. The case in which growing integrated negative energy is possible is when the 
squeeze parameter is small, r < 1. In this case, the energy density in an inertial frame is 
almost sinusoidal, with the positive energy outweighing the negative energy only in order r 2 . 
The effect of the periodic motion of the accelerated observer is to introduce Doppler shift 
factors which enhance the negative energy compared to the positive energy. The accelerated 
observer then sees the negative energy blueshifted and the positive energy redshifted. In 



the cases of perpendicular motion treated in Sect. II A and III A, the effect is a transverse 



Doppler shift, and is hence of order A 2 , where A is the oscillation amplitude. For the parallel 



motion cases in Sects. II B and III B , the effect is a linear Doppler shift, leading to an effect of 
order A. It is possible to have growing negative integrated energy density even for arbitrarily 
slow motion, which means arbitrarily small A. However, for a given A, the squeeze state 



parameter r is constrained by relations such as Eqs. (28) and (41), which limit the rate of 
growth. Note that non-relativistic motion is not a requirement for growing negative energy, 
and the numerically integrated results depicted in Figs. [T] and [3] are for relativistic motion, 
but small squeeze parameter. 

We studied a model which gives an operational meaning to integrated negative energy 



density in the form of defocussing of bundle of worldlines. In Sect. IV, we analyzed the 



Raychaudhuri equation for the expansion along a bundle of accelerated worldlines. The 
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motivation for this study is that positive energy leads to attractive gravitational effects and 
hence focussing, so negative energy should do the opposite. This expectation was born out in 
our results. However, the situation is complicated by the need to include an acceleration term 
in the Raychaudhuri equation, and the effects of the fluctuating velocity of the accelerating 
charged particles in a fluctuating electromagnetic field. In the cases which we examined, the 
gravitational effects and the acceleration effects have the same functional form. 

The effect treated in this paper bears a superficial resemblance to the effect treated in 
Ref. [20], which is a linearly growing or decreasing mean squared velocity of a charged 
particle undergoing sinusoidal motion near a mirror. The latter effect can be interpreted 
as non-cancellation of anti-correlated quantum electric field fluctuations. A charge at rest 
in the Casimir vacuum produced by the mirror is subjected to field fluctuations which can 
give or take energy from the charge for a time consistent with the energy-time uncertainty 
principle, but this effect will be cancelled by a subsequent anti-correlated fluctuation. The 
sinusoidal motion upsets this cancellation, and allows the mean squared velocity to grow or 
decrease, depending upon the phase of the oscillation. The effect discussed in the present 
paper also involves linear growth, but does not have an obvious interpretation in terms of 
non-cancelling fluctuations. The natural interpretation seems to be in terms of Doppler 
shifts which can be arranged to enhance negative energy and suppress positive energy. A 
topic for future research is to study further the connection between these two effects. 

Another topic is to understand to relation between the growing integrated negative energy 
and the general worldline quantum inequality of Fewster [TJ] . This inequality is difficult to 
evaluate explicitly for the sinusoidal worldline considered here. In this case, the inequality 
must be weak enough to allow the linear growth found here, but it might provide insight 
into the allowed behavior in situations more general than we have treated. 

A further question of interest is the possible physical consequences of accumulating neg- 



ative energy beyond those discussed in Sect. IV. A possible detection model for negative 
energy was proposed in Ref. [T7j, in which negative energy can suppress the decay rate of 
atoms in excited states. The atoms in this model are moving along inertial worldlines, but 
it might be possible to devise a more general model involving non-inertial motion. 

Let us also stress that our results do not in any way invalidate or diminish the implications 
of the quantum inequality bounds for inertial observers. The strength of a quantum inequality 
bound may depend on the particular observer chosen, but the validity of the bound does 
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not. As an example, suppose one is using a quantum inequality, applied to the motion 
of a particular inertial observer, to determine constraints on the geometry of a traversable 
wormhole. Let us further assume that in this case, the quantum inequality provides a very 
strong constraint. Now suppose one looks at the same problem from the point of view 
of, say, a different inertial or an accelerating observer and finds a much weaker bound. 
The weakness of the latter bound does not invalidate the strength of the previous bound. 
The observer whose motion provides the strongest quantum inequality bound implies the 
strongest constraint on the geometry of the wormhole. The latter cases simply yield true 
but weaker bounds. 
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